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Abstract 



We study neutron stars in a varying speed of light (VSL) theory of gravity in which the 
local speed of light depends upon the value of a scalar field </>. We find that the masses and 
radii of the stars are strongly dependent on the strength of the coupling between 4> and the 
matter field and that for certain choices of coupling parameters, the maximum neutron star 
mass can be arbitrarily small. We also discuss the phenomenon of cosmological evolution of 
VSL stars (analogous to the gravitational evolution in scalar-tensor theories) and we derive 
a relation showing how the fractional change in the energy of a star is related to the change 
in the cosmological value of the scalar field. 
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PACS numbers: 04.40.Dg, 04.50.+h 



1 Introduction 

Recently, there has been interest in the possibility that the speed of light c might have been 
larger in the past. The primary theoretical reason for considering this possibility is that such 
a variation in c could solve the horizon and flatness problems of big-bang cosmology, without 
needing to postulate the existence of an inflationary epoch in the early history of the Universe 
PP. This purely theoretical work was given added impetus more recently after the discovery by 
Webb et al. 2 that, according to observations of wavelength shifts in the absorption lines of 
distant quasars, the fine structure constant a := e 2 /{he) seems to have been smaller in the past. 
This is consistent with the assumption that c was larger in the past and, although more recent 
evidence [3] has constrained a variations more strongly than the Webb data, the possibility of 
a time varying a remains. 

One can assume that any changes in the value of a are due to variations in one or both of 
h and e. However, the physical consequences of allowing either of these to vary are different 
from those that arise from a varying c theory. For example, Avelino & Martins [3] have shown 
that the horizon and flatness problems can only be solved by assuming that c is varying. If 
instead one assumes that either of the other two constants is responsible for the variation of a, 
the horizon and flatness problems remain. 
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To analyse the consequences of allowing a physical "constant" to vary, one needs a self 
consistent theory which is capable of describing this variation and its effect on other physical 
processes. Since a, and by implication c, may change significantly only over cosmological time 
and distance scales, it is necessary to develop a theory of gravity which embodies such behaviour. 
This causes immediate problems when trying to formulate a varying speed of light (VSL) theory 
since the constancy of c is a fundamental assumption of relativistic physics. Several attempts 
have been made to reconcile such a variation with a relativistic gravitational theory, including 
a bi-metric theory developed by Clayton & Moffat [3] and Magueijo's scalar field theory [SJ[7]. 
It is the latter that we shall focus on here. 

In Magueijo's theory, the timelike coordinate in a coordinate chart no longer represents time 
as measured by a physical process and is instead replaced by a fourth length coordinate £. When 
described in terms of this new coordinate, the structure and properties of a manifold are the 
same as in any relativistic theory of gravity and all of the rules of differential geometry remain 
unchanged. In particular, the manifold has a metric, a Riemann curvature tensor, and so on, 
all of whose components can be expressed in a (£, x l ) coordinate basis, where the x % are spatial 
coordinates. Time as measured by physical processes is not in general usable as a coordinate on 
the manifold. However, once a solution is found, physically measured time intervals can then be 
determined using the relation dt = c~ l dQ, where c is treated as a scalar field on the manifold. 
The fact that t is no longer necessarily a timelike coordinate is apparent when one considers the 
possibility that in a solution, surfaces of constant c are not necessarily spacelike, so that V a c 
will not be timelike. The VSL theory is superficially similar to Brans Dicke (BD) theory except 
that the scalar field, instead of determining the strength of the gravitational field, governs the 
rate of physical processes and, by implication, the rate at which clocks record time. 

In any particular VSL theory, one would naturally expect spatial as well as temporal vari- 
ations in c in a spacetime with an inhomogeneous matter field. This is indeed the case, and 
Magueijo has analysed vacuum, spherically symmetric solutions representing both the weak 
field produced far from a central source and VSL black holes jjj. In this paper, we analyse the 
structure of neutron stars in Magueijo's VSL theory. The purpose of this is twofold: first, to 
determine the degree to which a VSL neutron star differs from its counterpart in general relativ- 
ity (GR) or BD theory and second, to serve as a starting point for the theoretical treatment of 
binary pulsars in a VSL theory, with the ultimate aim of using binary pulsar data to constrain 
the parameters of the theory. 

The plan of the paper is as follows. In Section |21 we give the action for the theory and derive 
the field equations for spherically symmetric spacetimes. In Section |3] we discuss the vacuum, 
exterior part of the solutions. In Section 0] we give the equation of state (EOS) of the fluid 
interior and discuss the issue of particle conservation, while in Section [5] we find the expression 
which describes the total energy of the solutions. In SectionElwe present the results of numerical 
integrations of the equations of structure and in Section [7| we discuss the evolution of the mass 
of a star in a cosmological setting. Finally, in Section |H1 we present a few concluding remarks. 



2 Action and Field Equations 

In Magueijo's VSL theory, the local speed of light c is determined by the value of a scalar field 
(f) such that c = Cooe^, where Cqo is the speed of light measured far from all gravitating sources. 
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In units in which G = 1, the action is given by 1 

I = ^T I d 4 x^g~ (n - Kg^VactN^) + IQ^lJ - ±- [ d 3 xVh e e<t >&, (1) 

107T JV L V ' J 07T JQV 

where V is an arbitrary volume of integration, dV its closed boundary, and h ab and the 
induced metric and extrinsic curvature on this boundary. The surface term guarantees that 
the variational derivative of / vanishes on dV and is important when determining the mass of 
an asymptotically flat spacetime. The constants e, (3 and k are arbitrary. In |S] the condition 
e — (3 = 4 was imposed since the units of the matter and gravitational terms of the action ([T]) 
differ by a factor of c 4 . We shall assume this condition holds here. In addition, we shall assume 
that the "minimum coupling" assumption of 6 holds so that L m is independent of (j). 
The generalised Einstein equations which lead from this action are 

G ab =U + e 2 ) V a 0V 6 - {Ik + e 2 )g ab V c <p\7 C <P + e(V a V fe </> - 5ab V c V c 0) + ^~ e ^T ab , (2) 



where T ab is the matter energy-momentum tensor, while the scalar field satisfies the wave equa- 
tion 

V a V a (t> = -eV>V a + ^ - (eT - 2pL m ) . (3) 

The local matter conservation law for the theory is 

V a T ab = -f3 (T ab - g ab L m ) V o 0- (4) 

In the vacuum case, where L m = and T ab vanishes, the action and field equations are identical 
to BD theory with 

u, = 4, $ = e ^, (5) 

where u and $ are the BD coupling parameter and scalar field. 

We are interested here in asymptotically flat, static and spherically symmetric solutions to 
the VSL field equations. We use isotropic coordinates and write the line element as 

ds 2 = -B d( 2 + A (dr 2 + r 2 d6 2 + r 2 sin 2 9\ , (6) 

where A and B are functions of r, and £ is a timelike coordinate with the units of length. We 
assume that the matter is composed of a perfect fluid with pressure p, energy density p and a 
fluid 4- velocity U a , for which the energy-momentum tensor has the usual form 

T ab = (p + p)U a U b +pg ab . (7) 

In the coordinates we are using, U a = {B~ 1 / 2 , 0, 0, 0) and is orthogonal to the hypersurfaces of 
constant £. 

To integrate the resulting field equations, it is convenient to introduce an auxiliary function 

1 dB 1 dA 
Z -2B^F + 2Ad^- (8) 



In terms of these variables, the field equations reduce to the following set of coupled ODEs: 



„„ 3^l' 2 2A' Arr A'\ Ak$ 2 „ 2 {( 

A" = — + etf [ AZ - — ) Z- + 8irA 2 e {e - 

AA r V 2 



p (e 2 - e_0) + 3pe 2 _ 
2k + 3e 2 9 



(9) 



The constants e and f3 we use here correspond to the constants a and b used by Magueijo |H||7j. We use the 
former notation to avoid confusion with tensor indices. 
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(10) 
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where the prime denotes 4-. 



3 Exterior Solution, Boundary Values and Scaling Relations 

For a given set of values of e, (3 and k, the solutions are parameterised by the central density 
Pq. As in GR, the surface of a stellar solution is defined to be where the pressure vanishes and 
at this point the solution is matched to a vacuum exterior. Since the vacuum field equations for 
this theory are identical to BD theory, the vacuum solutions are also identical and, in terms of 
the variables we are using here, are given by 

/ /v\ 4 tr - m\2(i-c/a-i/a) 

•^H) (t4) ' (13> 

Z = -4^—2 ( 14 > 



and 

ch = fli„ 4- 

eA \r + p, 

where the subscript "oo" denotes the value of quantities at spacelike infinity. In these equations 
p, C and A are constants and satisfy the constraint 

A 2 = (C + l) 2 -c(l-g). (16) 

Equations (|13|) and (|14|) together give the more familiar metric component 



B-D^[—^)'\ (17) 
\r + pj 



These equations hold only when e / 0. 

We shall assume that the interior neutron star solution matches smoothly to the above 
vacuum solution at the boundary of the star, in the sense that A, Z, (ft and the derivatives of 
these functions are continuous at the boundary. We can then find the values of the parameters 
p, C and A as follows. Differentiating eqns (|13ft to (|15[). one can show that 



X:=Z + e<f/= , 2//2 9 ■ (18) 
r° — p A r 



Rearranging this equation, we have 



Xr 3 

" = V4Tx? (19) 
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which allows us to determine the value of /x. From eqn (|14|). one can show that 



which allows us to calculate A once \i has been found. Finally, from eqn (|15[). we have 

C = ^!M, (21) 

which allows us to calculate C . 

The minimal VSL theory mentioned in [7j is the one with e = and its field equations are 
much simpler than in the more general theory. For this particular case, the above equations 
have to be modified. The functions A, B and Z are given by eqns (|13j) . (|14j) and (|17|) with the 
constant C = 0. The scalar field is now given by 

2 -« 2 \ 

(22) 




where g is a new constant, related to A and k by 

q 2 K = 2-^. (23) 

The Schwarzschild limit discussed in [7] is the special case of these equations in which A = 1. 
Note that in all cases in which e = 0, the <j) field need not be homogeneous. 

Magueijo has shown that, by carrying out a weak field analysis and interpreting the 
results within the framework of the PPN formalism, the parameters k, e and (3 are related by 
the equation 

7 - g2 + K + ^ (24) 

where 7 is the PPN parameter which governs, amongst other things, the deflection of light rays 
near to a massive body such as the Sun. Its limiting value is currently [S] 

|7- 1| < 0.0003. (25) 

Solving this eqn (|24j) for n we have 

e(e-(3) 



7-1 



+ \ef3 - 2e 2 , (26) 



which shows that k decreases as e — > for the case when e — (3 = A. For values of e of the order 
unity, eqn (|26|) restricts k to have a large value similar to that of uo in BD theory. However, if 
we allow e to have an arbitrarily small value, there is no minimum value to k. 



4 Equation of State and Particle Conservation 

We assume that the perfect fluid matter of the star consists solely of neutrons, each with rest 
mass m and local number density n. Denoting the internal energy per particle by II, the matter 
Lagrangian is then given by 

L m = -mc 2 n(l + n) = -p, (27) 
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the latter equality denning the relationship between n and p. Under the minimum coupling 
assumption outlined in 0, the particle rest energy mc 2 is constant and hence L m is independent 
oi</>. 

We shall consider two equations of state. The first is that of a a non-interacting gas of 
Fermions, which may be given in the parametric form 9 



V 



4 5 

m c 
24vr 2 /i 3 



X (2 X 2 ~ 3) vV + 1 + 3 log( X + vV + 1) 



(28) 



A c; 

m C 3 / 9 , 1 



This means that, in practice, eqn ()12() becomes a differential equation for \- Note that, under 
the minimum coupling assumption of jHj, both Ti and the product mc scale with in the same 
way, so that mc/U is independent of <p and hence so are numerical coefficients of eqns (128(1 . 

More realistic equations of state are generally given in tabulated form and require complex 
calculations which take into account particle phase transitions at various densities. However, 
polytropic equations of state fitted to the tabulated data do exist and our second EOS is of this 
type. The pressure and density are given in the parameterised form 



p = Kno fn[ — ) 



p = nm + 



n / 

Kno m ( n \ r 



(29) 



r-i \n 



where the constants are chosen to have the values no = 1.0 x 10 44 m -3 , K = 0.0195 and T = 2.34. 
The equations are parameterised by the particle number density n. This is the EOS used by 
Damour and Esposito-Farese in their investigation of strong field scalar effects in neutron stars 
jlUj and is based on tabulated data given in As before, we assume that the minimum 

coupling assumption holds so these parameters are independent of the local value of c. 

The first of our equations of state is far from realistic. However, we are interested here 
primarily in the effect a VSL theory has on the properties of a neutron star and we expect the 
main features of the solutions to be independent of the choice of EOS. In addition, our choice 
makes the results we find here easier to compare with the GR stars studied by Oppenheimer & 
Volkoff |12j and the BD solutions of Salmona The second EOS, being more realistic, will 
be used to place tentative constraints on the theory. 

Despite the presence of the factor of in the matter part of the action If!]), the conserved 
particle number N for a neutron star solution is the same as in a metric theory, namely 

N = f 8 4Trr 2 A 3 / 2 ndr, (30) 



o 



where n is the local particle number density. In the second EOS outlined above, this is the same 
quantity which parameterises those equations while in the simpler, non-interacting Fermion 
EOS, n and % are related by 



m°c 



n =(i^?J^ (31> 

In this relation, the bracketed combination of parameters is independent of c due to the minimal 
coupling assumption. 
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We note here that for a non-perfect fluid matter source, the conserved charge should in 
general include explicit reference to the fact that L m is not universally coupled. One can show 
that, for matter Lagrangians L m (tp) dependent on one or more matter fields tp which are invariant 
under a field transformation t/j — > t/j + 6i(), the Noether current associated with the symmetry is 
given by 

where I a is the corresponding current in a metric theory in which L m is universally coupled 
to the metric. This means that one should include a factor of in the expression for the 
conserved charge. Such a factor appearing in the Noether current would indicate a violation of 
a local energy conservation law. One can show that this does not occur for a perfect fluid source 
as follows. Dropping the assumption that the solutions are static and contracting eqn (j3J with 
Ub we have 



U b V a T ab = -pU b [(p + p)U a U b + (p + p)g ab \ V a <f> (33) 

which vanishes. Hence, at least for a perfect fluid, the non-spatial components of eqn Q are 
trivially zero and the right hand term contributes only to the pressure gradient. Thus there is 
no violation of the local energy conservation law in this case. 



5 Mass 

Since the solutions are spherically symmetric and asymptotically flat, one may write the ADM 
mass as the integral 

POO 

Madm = / ^r 2 G ab U a U b dr. (34) 
Jo 

This quantity is independent of the theory of gravity being considered and its value is determined 
solely by the metric. For both BD and VSL theories, from eqn (|13|) . the ADM mass is also given 

m adm = ?£+i>e. (35 ) 

The Poisson mass, discussed in (7j, may be defined by the relation 

M P := lim (\r 2 B'\ = ^, (36) 



-co 



A 



where the second equality follows from using eqn (|17|), 

In theories of gravity other that GR, Madm does not describe the total, physical energy of 
an asymptotically flat solution. For example, in BD theory the energy of a spherically symmetric 
solution is the Tensor mass |14j 



oo 



Madm - \ Hm r 2 — 



1 'M A dm-\Qs) (37) 



so called because it is the active gravitational mass which affects the motion of a test particle 
composed entirely of gravitational field self energy (such as a black hole). The second equality 
in eqn (|37|) defines the scalar charge Qs associated with <3?, which characterises the extent to 
which the gravitational energy of the star differs from that of a similar star in GR, and we use a 
tilde to denote quantities relating to BD theory. Equation (|3Tj) includes an overall factor of <&^~ 
to account for the fact that, as $ varies, the strength of the gravitational coupling G varies and 
the gravitational energy, which includes implicitly a factor of G, should vary in the same way. 
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To find an expression for the energy of a solution of the VSL theory, we shall find an 
expression for the Hamiltonian. The procedure is almost identical to that given for dilaton 
theories by Ho, Kim and Park ^S], so we only give the result here. We consider only a static 
spacetime; the generalisation to a time-dependent spacetime is as straight forwards as it is in GR. 
We assume that the spacetime is asymptotically flat and use the usual 3+1 splitting, denoting 
spacelike hypersurfaces by £ and the hypersurface orthogonal unit timelike vectors by U a . The 
timelike, 2-dimensional boundary of integration at spacelike infinity we denote by 5oo and its 
outward pointing, spacelike unit normal by n a . Starting with the action (^Q), decomposing 7Z 
into the Ricci scalar of the spacelike hypersurfaces and a divergence term, and integrating 
by parts several times, gives the VSL Hamiltonian 

H = - d 3 x^ (®R - (« + 2e 2 )h^d i( pd^ - 2eA(t>) + e^L r 

I d 2 x^B^e^ (0 - en a d a <t>) , (38) 

07TG JSoa 

where 9 is the extrinsic curvature scalar of, and a a b the induced metric on, the surface 5qo- The 
integral over £ vanishes by virtue of Hamilton's equation 5H/5\^B = 0, where \f~B is the lapse 
function, and we are left with the surface term. In general, this diverges and we must subtract 
a reference term H B . Following we choose H B to have the form 

H B = ^[ d 2 x^B^e t(t>B e B , (39) 
8vr J Soo 

where the subscript U B" denotes quantities evaluated when £ is intrinsically flat. The Hamil- 
tonian is then 

H=— f d 2 xVB^ e** (0 - en a d a (j)) - H B . (40) 



8vr J S , 

We assume that (f) is constant over Soo and denote its value by (j)^. We further assume the 
4>B = 0oo- One can then show that eqn lf4~0"j) reduces to 

H = (M ADM - \eQ s ) , (41) 

where 

Qs-=^- [ n a d a ^> (42) 
is the scalar charge associated with <j>. In a spherically symmetric spacetime, 

Q S = km (r 2d f) . (43) 

r^co \ dr J 

As one would expect from the form of the VSL action, eqn (|41|) is similar to eqn (|37jl. 

There is an ambiguity in the derivation of eqn (|41[) in that it is not immediately clear if one 
should multiply H by an additional power of e^°° to account for the effect on the total energy 
of varying the value of c throughout a solution. This is related to the problem discussed in [S] 
of deciding on the power of Cqo one must multiply the action Jy) by. To show that in fact eqn 
()41|) needs no such additional factor we consider the energy of a star in the weak field limit. 
Combining the second equality of eqn (|27j) with the definition of pressure p := mc 2 n 2 dll/ 'dn, 
one can show that 

(p + P) n ' = n p' ■ (44) 
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Combining this with eqn (J12[) and integrating gives 

mc 2 n = (p + p)y/B e w . (45) 

We can define a "rest energy" Mr = mc 2 N for the star, which is simply the sum of the energies 
of its constituent particles and is its total energy in the weak field limit. Combining eqns (|3UI) 
and (|4*5|) we have 

POO 

M R = 47rr 2 VB A^ 2 e^{p + p) dr. (46) 
Jo 

In the lowest order weak field limit, p » p, g a b — ► r] a b and <p is constant. In this case, eqns ([%]) 
and ()34jl together imply that the weak field ADM mass is 

roo 

Madm = / 4^r 2 / oe (/3 " e)</, °° dr, (47) 
J o 

while the Hamiltonian in this limit is simply e^ 00 Madm- Under the lowest order weak field 
approximation, this is identical to the quantity defined by eqn Q46|) indicating that no additional 
powers of e^°° need be included in the energy definition and so we shall adopt eqn 1)4 ljl as the 
definition of the energy of the star. Due to this expression's similarity with eqn 1)37)1. we shall 
henceforth denote the energy of VSL solutions by Mt- 

We may check that our definitions of energy and particle number, defined by eqn ()3())) are 
consistent as follows. A necessary condition for the stability of stellar solutions in any theory of 
gravity is that, given a set of solutions parameterised by a smoothly varying quantity such as 
their central density, the extremal values of the total energy and the conserved particle number 
should coincide Using a method similar to the one outlined in JH] for boson-fermion stars, 
it is relatively easy to prove that N defined by eqn ()3U|) and our choice of eqn 1)41)1 for the energy 
of a solution are consistent with this requirement in that these two quantities do have coinciding 
extrema. 

In a metric theory, Mp defined by eqn 1)36)1 is also the Kepler mass Mk of the star. However, 
this is not the case here: the Kepler mass in a VSL theory includes a term proportional to (3 
which embodies the non-minimal coupling between L m and the metric and, in terms of the 
variables we are using here, is given by [Jj 

M K = M P + PQ S , (48) 

where, in terms of Madm and Qg, the Poisson mass is 

M P = Madm - eQg. (49) 

These definitions are similar to those in BD theory except that in the later case, Mp and Mk 
are identical. 

6 Numerical Results 

Equations © to ()12)) are invariant under the rescaling 

+ A -»■ Ae^-ft". (50) 

We use this property to set (j)^ = for all solutions. In addition, eqns © to 1)12)) are invariant 
under the rescaling 

A -» a 2 A, r-*-, Z -> aZ (51) 



9 



and we use this to set A^ = 1 for all solutions. Note that this last rescaling implies that 
fi — > fi/a. The function Z automatically satisfies the boundary condition = 0. We calculate 
fi, C and A for each solution using eqns (|19|) to (|21|) . in which we evaluate A, Z, <ft and their 
derivatives at the outermost point of the interior solution. Equation ((Td]) is then used to check 
the consistency of the numerical calculations. We characterise the size of each solution by the 
Schwarzschild radius R, which is related to r by 

R:=^JT s r s , (52) 

where the subscript s denotes quantities evaluated at the surface of the star. This is the same 
definition of radius used in |13j . 

Figure 1 shows curves of mass against central density for several sets of solutions with EOS 
1281 For each set, the value of e is fixed, we impose the condition e — (3 = 4 and choose the 
minimum value of k consistent with eqn (|26|) . The mass of each solution is calculated using eqn 
(|41|) . The e = 4 solutions are indistinguishable from those in GR since they are identical to 
solutions in BD theory with a large value of u. The Figure shows that as e and k decrease, the 
energy of all of the solutions decreases. Masses Mt and Mk, along with other data from the 
same sets of solutions, are shown in Table we do not show values of Mp since they are very 
close to corresponding values of My. The data in the table indicate that, as well as the mass, 
the radii of the solutions also decreases rapidly with e. 



e 


K 


K/e 2 


M T 


M K 


Qs 


R 


4.0 


53290 


3330 


0.711 


0.712 


-0.0000 


10.1 


1.0 


13330 


13330 


0.699 


0.712 


-0.0036 


10.0 


0.1 


1333 


1.33 x 10 5 


0.588 


0.711 


-0.0311 


9.50 


0.005 


666 


2.67 x 10 5 


0.495 


0.704 


-0.0524 


8.99 


0.001 


13.3 


1.33 x 10 7 


0.192 


0.598 


-0.1017 


6.62 


0.0005 


6.66 


2.67 x 10 7 


0.096 


0.502 


-0.1016 


5.27 


0.0001 


1.33 


1.33 x 10 8 


0.012 


0.269 


-0.0642 


2.66 



Table 1: VSL neutron stars with e — f3 = 4 and central density po = 2.9 x 10 18 kgm~ 3 . The 
energy Mt is calculated using eqn (flTj) , the Kepler mass Mk using eqn 1)48(1 and Qs is the scalar 
charge associated with <p. All three are measured in units of solar mass. The radius, in units of 
km, is calculated using eqn (|52[) . 

The response of the properties of the VSL stars to changes in the scalar field coupling strength 
is in marked contrast to to that for neutron stars in BD theories, data for which are shown in 
Tabled (Further data on the structure of BD stars can be found in Salmona |13j but note that 
the mass used in this reference is the Poisson mass and the boundary value of <3? is different from 
our value of = 1.) Values of u have been chosen in the Table to match those of k in Table 
Hand we show masses Mt and Mk', the Poisson mass Mp is identical to Mk- The data show 
that increasing the strength of the scalar field's coupling has a moderate effect in decreasing the 
energy and radius of the stars, while increasing the Poisson mass. 

There are two reasons for the differences between BD and VSL neutron stars' properties. 
From eqn (|11|). for small values of e, the strength of the coupling between </> and the matter 
varies as («/e) . This quantity increases as k decreases, which means that the (j) field becomes 
more inhomogeneous and the magnitude of Qs increases as we decrease e. This behaviour is 
very similar to that for BD theory. However, in contrast with BD theory, the ADM mass, 
Poisson mass and energy Mt differ by multiples of eQs and the product of these two quantities 
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M T 


M P 


Qs 


R 


1333 


0.711 


0.711 


-0.0004 


10.1 


666 


0.711 


0.712 


-0.0008 


10.1 


13.3 


0.705 


0.723 


-0.0365 


10.1 


6.66 


0.699 


0.732 


-0.0655 


10.0 


1.33 


0.673 


0.764 


-0.1818 


9.95 



Table 2: Neutron stars on BD theory with central density po = 2.9 x 10 18 kgm -3 . The energy 
Mt is calculated using eqn (|3*Tj) , the Poisson mass Mp using eqn (|3*B^) and Qs is the scalar charge 
associated with the BD scalar field All three are measured in units of solar mass. The radius, 
in units of km, is calculated using eqn (|52j) . For BD theory, the Poisson and Kepler masses are 
identical. 



decreases as e decreases. This is in contrast with BD theory, where the larger value of Qs serves 
to increase the Poisson (or Kepler) mass to a value significantly larger that other masses. 

The second, and more significant, cause of structure difference between the two theories is 
due to the non-universal coupling between the matter and the metric in the VSL action. From 
eqn (|12ft . the pressure gradient has a term involving the product /3cf>' . For all of these solutions, 
<fi' < and, under the assumption that e — (3 = 4, (3 < 0. Hence the product f3<f>' > and since 
p must be a decreasing function of r, eqn 1)12(1 implies that the term involving (3 increases the 
magnitude of the pressure gradient, leading to p — * at a smaller value of R than for a star 
in metric theory with an identical central density. The significantly smaller physical size of the 
star means that its energy is also much smaller. 

Figure 2 shows mass against central density for sets of solutions with EOS 1291 with the 
same limits on the parameters e and (3 as for the solutions of Figure 1. These solutions show 
the same overall behaviour as those with the simpler EOS: as the parameter r) decreases, the 
solutions become less massive for a given central density. As before, the e = 4 solutions are 
indistinguishable from those of GR. In addition, the central density at which the maximal mass 
solution occurs increases as e decreases, again in a similar way to the solution shown in Figure 
1. In addition we have found that, for a given central density, the radius of a solution decreases 
as e decreases. However, both this feature and the reduction in mass are less severe than for 
the simpler EOS and both maximal mass and minimum radius of the maximal mass solution 
decrease far less rapidly with decreasing e. This is due to the difference in the stiffness between 
the two equations of state. 

Since the Figure 2 solutions have a reasonably realistic EOS, we can use them to place a 
limit on the parameters of the VSL theory for the case studied here in which e — j3 = 4. The 
masses of neutron stars in binary systems detected to date have masses in the range (1.36 ±0.8) 
solar masses, while neutron stars in orbit about other types of companion body have masses 
of at least 1.4 solar masses (for a review of this data and references, see [2^). Hence, for a 
VSL theory to be viable, it must allow for the existence of neutron stars whose mass is at least 
1.3 solar masses. From Figure 2, this implies that the parameter e must satisfy the inequality 
e > 0.00015. 
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7 Sensitivity and Evolution of Mass 



For many theories of gravity, such as BD theory and other more general scalar-tensor (ST) 
theories, solar system observations place strong constraints on the theories' parameters via the 
PPN formalism, while limits derived from observation of strong field phenomena are far weaker. 
However, for VSL theories with small values of e, the PPN limit (|26|) on k is extremely weak 
and we should therefore look to strong field observations, in particular binary pulsar data, to 
provide the strongest constraints. The situation here is similar to that found by Damour and 
Esposito-Farese for certain classes of ST theory which satisfy PPN limits to an arbitrarily high 
degree but show significant strong field effects |lf)j and references therein. We shall leave the 
full treatment of the motion of binary pulsars in a VSL theory to a future publication. Here 
we shall briefly mention some of the key ideas involved and use part of the formalism to discuss 
how VSL stars will evolve in a cosmological setting. 

The analysis of binary pulsar motion is based on the parameterised post-Keplerian (PK) 
formalism, which treats each star as a point mass moving according to quasi-Newtonian equa- 
tions of motion. Details of the particular theory under study are made manifest by allowing 
the active, passive and inertial masses of each star to vary according to the value of any am- 
bient, non-metric gravitational fields. At the lowest level of approximation, the variability of a 
particular mass is quantified by the sensitivity s, given by 



y, d{\OgM) 



<9(log tjj A ) 



(53) 

■ oc 



where M is the mass under consideration and the ip are the non-metric gravitational fields. 
In ST theories of gravity, there is only one extra field, the scalar field and it turns out that 
active, passive and inertial masses are proportional to each other. For the VSL theory we are 
considering here, we expect much of the formalism to be very similar to that for BD gravity, 
since the vacuum, inter-body field equations are the same for both theories. Since the BD and 
VSL scalar fields are related by $ = e^, we shall define s to be 

_ g(log M) J_dM_ 
S ~ d{e4>) ~ eM d(j) ' ( ' 

where again M denotes any mass. 

We shall next derive an expression for the sensitivity associated with the Hamiltonian, iden- 
tical to the VSL tensor mass Mf, of a compact object in a VSL theory. We consider first the 
variation of the Hamiltonian (|38j) induced by a variation of (j) on Soo. Taking the functional 
derivative of this equation, performing several integration by parts and noting that the volume 
term vanishes and the terms involving d a (5(p) in the surface integral cancel, we are left with 

5H = -^— [ d 2 xVB^e^ (eO - en a U b V b U a + nn a d a A 5^ + 5H B , (55) 
ovrG J Soc v ' 

where we have assumed that 5<j) is constant over and have denoted it by <5(/>oo- The reference 
term Hb is linear in e^°° . Thus 

5H B = eHgS^oo. (56) 
Combining the above two equations with eqn (|38|). we have 

5H = -5- / d 2 xVB^e^ Uk + e 2 )n a d a <j) - en a U b V b U a ] <50oo + eHd^. (57) 

07TG J Soo 1 J 
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Using the definition ()36|). one can show that 



tfxyfBo n a U b V b U a = \M P . (58) 



8vrG J Sa 

In addition, the first term of eqn ()57|) is simply ^(k + e 2 )Qs- Thus we have 



5H = 



UK + e 2 )Qs-\eMp\5ct> 00 + eH5<t )oo . (59) 



Finally, using the relations (|41[) and (|49|). and taking the limit <5^>oo ~~ ► 0, one can show that 

= + e *- ____ ] g 5 . (60) 



4 

Hence the sensitivity associated with both H and the identical quantity My is 

Due to the similarities between the two theories, we expect s to appear in the PK formalism for 
a VSL theory in much the same way as in a ST theory. Thus we expect that, for any particular 
mass M, presence of the 4> field introduces a strong field correction of the form M — > M(l + i/s), 
where v is a parameter which depends the type of mass considered. Equation ()61l) then indicates 
how strong field gravitational effects depend upon the values of e and /3. 

Equation 1)6 1|) may be used to investigate the gravitational evolution of a compact body in 
a VSL theory. The scalar field far from a compact body such as a neutron star should match 
smoothly to the cosmological value of the field. Hence, assuming that a body embedded in a 
cosmological solution can be modelled approximately by an isolated, asymptotically flat solution, 
one expects that the boundary field (j)^ to follow the cosmological evolution of <p. Equation 1)61)) 
then implies that the evolution of the mass of the body depends strongly on the relative values 
of e and k. For all of the solutions shown in Figure ^ and Tabled except the e = 4 solutions, 
s < 0. The primary reason for this is that, as e decreases, the ratio k/e which dominates 
the second term of eqn ()61j) increases dramatically. In a cosmological solution, <p should be a 
decreasing function of cosmological time so as to be consistent with the requirement that c was 
larger in the past. Use of eqn ()61|) then implies that the energy of these VSL stars increases with 
cosmological time, the largest rates of increase coinciding with smallest values of e. This is in 
marked contrast with ST theories for which, in general, the energy of a compact body decreases 
with cosmological time ^1 120] . 



8 Conclusions 

We have examined the properties of neutron stars in one particular VSL theory in which the 
local value of the speed of light is determined by a scalar field <p. The theory bears a more than 
superficial resemblance to BD and other ST theories, and we have found that the expression 
giving the energy of an asymptotically flat solution to the VSL field equations is very similar 
to it counterpart in ST theories. However, despite these similarities, the properties of VSL 
neutron stars are very different to their counterparts in VSL theories. In particular, for certain 
choices of e and 0, VSL neutron stars are limited to being much smaller those in GR and the 
dependence of mass on the strength of the coupling between <j> and the matter is extremely 
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strong. This allows one to use the mere existence of neutron stars to place constraints on the 
theory's parameters, even without having to analyse the motion of binary pulsars: the maximum 
neutron star mass allowed by the theory must obviously be at least as large as those of observed 
neutron stars. Unfortunately, neutron star masses are fairly sensitive to the EOS used (see, for 
example, Y22\). although we have placed a rather crude limit on the parameter e using a 
reasonably realistic equation of state. 

We have also given a brief analysis of the evolution of a VSL star in a cosmological setting 
found that, as the cosmological background evolves and 4> decreases in value, the energy of a 
VSL star increases in general. This result is independent of the EOS used or even of the type 
of matter of which the star is composed and could have serious implications for the stability of 
a star since its fractional binding energy (My — N)/N could eventually become positive. 

To properly place limits on the parameters of a VSL theory, one needs a full analysis of the 
predicted behaviour of binary pulsar motion. However, this is not the only area of astrophysics 
in which a VSL theory might give significantly different predictions to those of other theories. 
For example, it was shown in [7j that black holes in this theory are very different to their GR 
counterparts in that the proper time taken for an in-falling particle to reach the horizon can be 
infinite. If this is indeed the case, a varying speed of light may have a significant effect on the 
physics of black hole accretion disks since matter orbiting the black hole would never enter the 
black hole horizon. 
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Figure 1: Mass curves for neutron stars in the VSL theory with the simple, non- interacting 
Fermion EOS (|28|). The curves are labelled by their value of e. For each curve, e — (5 = 4 and 
the value of k is the minimum value consistent with the weak field, solar system constraints. 
The central density is given in units of kgm~ 3 . 
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Figure 2: Mass curves for neutron stars in the VSL theory with the realistic EOS ()29|) . Units, 
labels and curve parameterisation is the same as in Figure ^ 
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